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Abstract 



The model of induced quark currents formulated in our recent paper (Phys. 
Rev. D51, 174) is developed. The model being a kind of nonlocal extension of the 
\ bosonization procedure is based on the hypothesis that the QCD vacuum is realized 

■ by the (anti-)self-dual homogeneous gluon field. This vacuum field provides the 
analytical quark confinement. It is shown that a particular form of nonlocality of 

\ the quark and gluon propagators determined by the vacuum field, an interaction 

of quark spin with the vacuum gluon field and a localization of meson field at the 

■ center of masses of two quarks can explain the distinctive features of meson spec- 
trum: Regge trajectories of radial and orbital excitations, mass splitting between 
pseudoscalar and vector mesons, the asymptotic mass formulas in the heavy quark 

• limit: Mqq — > 2mQ for quarkonia and Mq^ mq for heavy-light mesons. With 

"fH , a minimal set of parameters (quark masses, vacuum field strength and the quark- 

CL|| gluon coupling constant) the model describes to within ten percent inaccuracy the 

Q^' masses and weak decay constants of mesons from all qualitatively different regions 

O . of the spectrum. 

^ PACS number(s): 12.39.-x, ll.lO.Lm, 12.38.Aw, 14.40.Gx 



H 1 Introduction 



Achievements of the Nambu-Jona-Lasinio (NJL) model in description of meson masses, 
decay constants and so on are well-known [|I|, |^, ^. This success can be explained by the 
bosonization procedure which makes possible to extract collective modes and dynamical 
breaking of SUl(3)xSUr(3) and Ua(1) symmetries. At the same time, an incorporation 
of quark confinement into consideration, description of heavy quarkonia and heavy-light 
mesons, radial and angular excitations of mesons, as well as different form-factors require 
an essential modification of the NJL model (e.g., see [^]). Another general disadvantage 
is the nonrenormalizability of the local four-fermion interaction. 

In recent paper we have suggested a model that in some sense can be considered as 
an extension of the standard NJL model. There are two crucial modifications. First, our 
model is based on the hypothesis that the QCD vacuum to be realized by the (anti-)self- 
dual homogeneous background gluon field. Second, the effective quark-quark coupling 
is described by the nonlocal four-quark interaction induced by the one-gluon exchange 
in presence of the (anti-) self-dual homogeneous gluon vacuum field. This vacuum field 
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ensures the analytical quark confinement and breaks the chiral symmetry. The model of 
induced quark currents gives a basis for investigating of all the above-mentioned problems 
/^from a general point of view. The main features of the model are as follows 0. 

- There is the quark confinement. The quark propagator being an entire analytical 
function in the complex momentum plane has the standard local ultraviolet behavior in 
the Euclidean region, and is modified essentially in the physical, i.e., Minkowski region. In 
contrast to purely chromomagnetic and chromoelectric configurations, the (anti-)self-dual 
homogeneous field is a stable configuration [0. 

- The one-gluon exchange is decomposed into an infinite sum of current-current in- 
teraction terms, in which the quark currents are nonlocal, colorless and carry a complete 
set of quantum numbers including the orbital and radial ones. This effective quark-quark 
interaction generates a superrenormalizable perturbation expansion. 

- The bosonization of the nonlocal four-quark interaction leads to ultraviolet finite 
effective meson theory. Mesons are treated as extended nonlocal objects. 

- The model contains the minimal number of parameters: the quark masses, quark- 
gluon coupling constant and the tension of the background gluon field. 

It was shown also, that the spectrum of the radial and orbital excitations is equidistant 
for sufficiently large angular momentum £ or radial quantum number n. In the heavy 
quark limit the mass of quarkonium tends to be equal to sum of the masses of constituent 
quarks. 

In paper ^ the main attention was paid to mathematical details of obtaining nonlocal 
quark currents induced by the one-gluon exchange in presence of the vacuum field and to 
formulation of the bosonization procedure based on these nonlocal currents. A motiva- 
tion and connection of the model with QCD has been discussed as far as possible. The 
present paper is concentrated on further development of the model and on application to 
systematic calculations of the weak decay constants and masses of mesons from different 
regions of the spectrum: the light mesons and their excited states, heavy quarkonia and 
heavy-light mesons. 

With the minimal set of parameters (quark masses, vacuum field strength and one 
quark-gluon coupling constant) the model describes all qualitatively different regions of 
meson spectrum to within ten percent inaccuracy. The reasons driving this successful 
description can be easily recovered. 

An interaction of quark spin with the vacuum gluon field, contained in the quark 
propagator, breaks the chiral symmetry and gives rise to the splitting between masses of 
the pseudoscalar and vector mesons with identical quark structure (p — vr, K — K*). This 
spin-field interaction drives also the weak decay of pion and kaon. 

Furthermore, the vacuum field produces three rigid asymptotic regimes for the spec- 
trum of collective modes. The spectra of radial and orbital excitations of light mesons are 
equidistant for £ ^ 1 or n ^ 1, i.e., they have Regge character. This is due to the specific 
form of nonlocality of the quark and gluon propagators determined by confining properties 
of the vacuum gluon field. After all, one concludes that the confinement is responsible for 
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Regge trajectories. Localization of meson field at the center of masses of a quark system 
provides other two asymptotic regimes. In the limit of infinitely heavy quark, a mass 
of quarkonium tends to be equal to a sum of the masses of constituent quarks, while a 
mass of heavy-light meson approaches the mass of a heavy quark: Mqq 2mQ — Aqq, 
Mqq — > rriQ + Aq^. The next-to-leading terms Aqq and Aqg do not depend on the heavy 
quark mass. The same reasons provide the correct asymptotic behavior of the weak decay 
constant for the heavy-light pseudoscalar mesons: /p ~ 1/ ^Jtuq. 

One can conclude that the (anti-) self-dual homogeneous background gluon field deter- 
mines rather definitely a behavior of the masses and weak decay constants in all different 
regions of the meson spectrum, and this behavior is quantitatively consistent with exper- 
imental data. 

Technically, these results are based on a decomposition of the bilocal colorless quark 
currents into a series of nonlocal currents with complete set of quantum numbers: spin, 
isospin, radial and orbital numbers. Tensor structure of these nonlocal currents is rep- 
resented by the irreducible tensors of the four-dimensional Euclidean rotational group, 
while their radial part is determined by a specific form of gluon propagator in the exter- 
nal vacuum field. This decomposition provides a new point of view on the renormalization 
problem: the Feynman diagrams appearing in each order of perturbation theory are ul- 
traviolet finite due to the nonlocality of the meson-quark interaction. 

The paper is organized as follows. In Sect. II we review the main points of the model 
with some modifications, that relate to choosing the point of localization of meson field and 
description of super-fine structure of the spectrum. These modifications reflect possibili- 
ties missed in f^. In Sect. Ill we consider the masses of light mesons and their excitations, 
heavy quarkonia and heavy-light mesons, as well as the weak decay constants. The basic 
approximations of the model and problems for further investigations are discussed in the 
last section. 



2 The model of induced nonlocal quark currents 

2.1 Basic assumptions, approximations and notation 

The representation of the Euclidean generating functional for QCD, in which the gluon and 
ghost fields are integrated out, serves as a starting point for many models of hadronization. 
Dyakonov and Petrov obtained this representation for the case of nontrivial vacuum gluon 
field m (see also 0) 



Z = J do^^c J DqDqexp^^J d^x ^ g/(x)(z7^V;, - m/)g/(x) + ^ 

where Np is the number of flavors corresponding to the SU{Np) flavor group and 

9'^ 



(1) 



3 



f 

The function G^^ "^" is the exact (up to the quark loops) n-point gluon Green function in 
the external field B'^. We will investigate the mesonic ((3'g)-collective modes and consider 
Eq. (|l|) with the quark-quark interaction truncated up to the term L2 

Z = J dcTvac J DqDqexpi^J rf^x ^ g/(x)(z7^V^ - m/)g/(a;) 

+ y // d'xd% j;{x)G%{x,y \ B)jl{y)^. (2) 

Representations ([l|) and (H) imply, that there exists some vacuum (classical) gluon field 
-B°(x), which minimizes the effective action (or effective potential) of the Euclidean QCD. 
In the general case, the vacuum field depends on a set of parameters {(Xvac}, and the 
measure (io"vac averages all physical amplitudes over a subset of {cTvac}, in respect to which 
the vacuum state is degenerate (for more details see and references therein). 

The quark-gluon interaction both in Eqs. (0) and (^ are local, and a decomposition 
over degrees of generates a renormalizable perturbation theory. It means, that an 
appropriate regularization is implied in Eqs. (|l|) and (@). This point has to be stressed here, 
since our final technical aim is a transformation of the interaction term in Eq. (^, which 
generates completely new superrenormalizable perturbation expansion of the functional 
integral (0). 

Let us identify all ingredients of these general formulas for the particular case of an 
homogeneous ( ant i-) self- dual vacuum field 

B^{x) = nB^{x), B^{x) = B^^x^, 
n = n^t", = n^n" = 1. 

The constant tensor B^^, satisfies the conditions: 

B^u By^^ B i^ipB pu B 



B^u ^uafjBo^fj T^SlB^jj^ 



where i? is a gauge invariant tension of the vacuum field. Since the chromomagnetic H 
and chromoelectric E fields relate to each other like H= ±E, two spherical angles (</?, 9) 
define a direction of the field in the Euclidean space. In the diagonal representation of 
h = n°'t°', an additional angle ^ is needed to fix a direction of the field in the color space 

n = t^cos^ + t^sin^ , < ^ < 27r. 
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The one-loop calculations and some nonperturbative estimations of the effective po- 
tential for the homogeneous gluon field argue (do not prove) that the potential could 
have a minimum at nonzero value of the field tension B (e.g., see 0, |TU|). We 
will assume, that the field under consideration realizes a nonperturbative QCD vacuum, 
and study the manifestations of this field in the spectrum of collective modes. Since the 
effective potential is invariant under the Euclidean rotations, parity and gauge transfor- 
mations, this vacuum should be degenerate with respect to the directions of the field in 
the color and Euclidean space and should be the same for anti-self- and self-dual config- 
urations. According to this argumentation, the field 5° in (|^) corresponds to the tension 
B minimizing the effective potential, and the measure da^ac has the form 

^ =*= 

where the sign X]± denotes averaging over the self- and anti-self-dual configurations. To 
simplify calculations and to clarify the technical side of bosonization procedure in presence 
of the background field, we will omit the integral over ^ in (^ and fix a particular vector 
n°- = S""^. In the fundamental (matrix t^) and adjoint (matrix C®) representations of 
SUc{3) one gets 

n = t' = diag(-^, --^) , B^,B,, = -{^fBH^,, 

n = C'' = , B^^B^, = -^K'B\, , 

7^54 = -i^45 = Kre = ~K^j = i , = diag(0, 0, 0, 1, 1, 1, 1, 0). 

The rest of elements of the matrix K are equal to zero. It is convenient to define the mass 
scale = ^B: 

/I 1 2 

B^pBp^ = -v^A^6^^, V = diag ( 3 

We would like to stress, that the averaging over directions of the background field in the 
color space should be incorporated into the formalism, and its role should be analyzed. 
But first of all, we would like to go as far as possible with the formalism, that is simpler as 
possible, and test, how the technique proposed in works in the meson phenomenology. 

2.2 Quark and gluon propagators 

The quark propagator Sf{x,y \ B) in Eq. (0) satisfies the equation: 

{'i'lpSp - mf)Sf{x,y I B) = -5{x - y), 
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and can be written in the form 



Sf{x,y I B) = ei'-''^y^Hf{x - y \ B)ei^^^y\ 



HAz) 



Hfip I B) 



V^{z) + m} + (aB) 
1 



6{z) 



(4) 



2vA 



1 ^ 
dte 2ua2 



l + t 



X 



P+ + R 



1 + e 



^(7/7)^ 



(5) 



where 



^± = ^(1 ± 75), = ' (^^2/) = x^B^^y^, , 

The function Hf is the Fourier transformed Hf. The upper (lower) sign in the matrix P± 
corresponds to the self-dual (anti-self-dual) field. 

The term {a 13) = <7a/3Ba/3 in Eq. (|^) (the second line in Eq. (|)) describes an interaction 
of a quark spin with the background field. One can see, that this spin-field interaction 
leads to the singularity l/nif for mj — » 0, which is a manifestation of the zero mode 
(the lowest Landau level) of the massless Dirac equation in the external (anti-) self-dual 
homogeneous field. The mathematical point is that the spectrum of the operator jfj.dfj, 
is continuous, whereas the spectrum of the operator 'j^'V^{x) is discrete and the lowest 
eigen number is equal to zero. Simple calculations give for m/ 



Hf{p I B) = 2e 



V 



1 1 

— + — (Pm7m + ^(7/p) 
iTLf mj 



Pt - 7(7/7) 



+ 0(1), (6) 



and 



lim lim m/-(o/-(a;)of (x + £:)) r = — lim lim mfTTHfielB) 

= - [ -r-^ lim TUfTrHfip I B) = — ^rA^ (7) 

Due to the spin-field interaction the quark condensate is nonzero in the limit of vanishing 
quark mass. This indicates that the chiral symmetry is broken by the vacuum field in 
the limit mj — (see also 0). It will be clear below that just the spin-field interaction 
gives rise to the splitting between the masses of the pseudoscalar and vector mesons and 
provides a smallness of pion mass. 
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In terms of the variable ( = p^"^^ the propagator i//(C | B) is an entire analytical 
function in the complex (^-plane. There are no poles corresponding to the free quarks, 
which is treated as the confinement of quarks. The following asymptotic behavior takes 
place: 

\ 0(exp(^)) = 0(exp(^)) if C - ±.oo (p^ - -00). ^ ^ 

Equations (H) shows the standard local behavior of the fermion propagator in the Eu- 
clidean region [j? — >■ 00), while in the physical region [jp' — > —00) we see the exponential 
increase typical for nonlocal theories (for more details about the general theory of nonlo- 



cal interactions of quantized fields see [TT[|). Below, the absence of the poles and the 



exponential increase will be referred as the confinement properties of a propagator. 

Function D^|j,(a;, y\S) in representation is the exact gluon propagator for the pure 
gluodynamics in presence of the vacuum field 5^. This function is unknown, and some 
approximation has to be introduced. For instance, the local NJL model corresponds to the 
choice = b^-^b^ybix — y). We go beyond this approximation and replace the function 
by the confined part 



D^,{x,y I B) = 6^,K\^^-^y^D{x - y \ X^)ei^-^y\ (9) 



of the gluon propagator 

Gfy{x, y I B) = D%{x, y \ B) + y \ B), 

which is a solution of the equation (for details see [^) 

(V^5^^ + AiB^y)Gyp{x, y\B) = -5^p5{x - y). 

The Fourier transform of the function D{z \ h?) is an entire analytical function in the 
momentum space. It has the local behavior in the Euclidean region, but increases ex- 
ponentially in the physical region. This function describes a propagation of the confined 
modes of the gluon field. Other terms -R^t' ^^^^ contain a contribution of the zero modes 
and an anti-symmetric part, will be omitted. 

Thus, the central point of our extension of the NJL-model consists in taking into 
account the confining influence of the background field both on the quark and gluon 
propagators. 

2.3 Color singlet bilocal quark currents 

Substituting gluon propagator (|TUp to the interaction term in representation (^, using the 
Fierz transformation of the color, flavor and Dirac matrices, and keeping only the scalar 
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J"*^, pseudoscalar J°^, vector J"^ and axial- vector J°"^ colorless currents, we arrive at 
the expression ||^ 

^2 = yECj// d'xdSj'\x,y)D{x-y\ K')J'\y,x), (10) 

bj 

J'^ix, y) = qf{x)M'ff,r'e^^^^ykf{y) , (11) 

r^ = i, r^ = ^75, r^ = 7^, 1^ = 757/., 
1 1 

= Cp = - , = = — . 

Here Mjj-, are the flavor mixing matrices {b = 0, A^^^. — 1) corresponding to the SU(A'^p) 
flavor group. In the case of SU(2) and SU(3) they are given by the matrices and 
respectively. 

Due to the phase factor exp{i{xBy)}, bifocal quark currents (0) are the scalars under 
the local gauge transformations 

q{x) ^ e-^^(^')g(x) , q{x) g(x)e^^(^) , 

Bu e-^'^(^)5„e*'^(^) + le-^'^(^)9.e^'^(^\ (12) 

9 

Let us transform integration variables x and y in Eq. ([T0|) to the coordinate system 
corresponding to the center of masses of quarks qf{x) and qf'{y) 

X ^x + ^fy, y ^x- if>y, = ^ , = — ^ . (13) 

Corresponding transformation of the quark currents looks as 

J'^ix^y) = qf{x)M)f,T'e'^^^y\f.{y) ^ 
qf{x + ify)M)^,T'e^^^y\f.{x - ^f>y) = 

qf{x)M'}f,r-^e^y^ff''^''kf'{x) = J^\x,y), (14) 
where V//' is a linear combination of the left and right covariant derivatives 

V//' {x) = ^f(d +^B{x)) - ^jid -iB{x)). 

These covariant derivatives indicate, that the currents ([1^) are nonlocal and colorless. 
Interaction term (|l^) takes the form 

= I^E^^ // rfWy^expj-^l J^^(x,y)/-^(y,x), (15) 

where we have made use of the representation (|^). The currents are deflned by Eq. ([T^). 
Transformation (|T3p turns out to be crucial for simultaneous description of the light-light, 
heavy-light and heavy-heavy mesons. 
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2.4 Decomposition of bilocal currents 

An idea of our next step consists in a decomposition of the bilocal currents (0) over some 
complete set of orthonormalized polynomials in such a way, that the relative coordinate 
of two quarks y in Eq. (^) would be integrated out. One can see, that a particular form 
of this set is determined by the form of the gluon propagator (exp{— in Eq. (|15])). 
The propagator plays the role of a weight function in the orthogonality condition. The 
physical meaning of the decomposition consists in classifying a relative motion of two 
quarks in the bilocal currents over a set of radial n and angular I quantum numbers. In 
other words, according to general principles of quantum mechanics the bilocal currents 
have to be represented as a set of quark currents with definite radial n and angular £ 
quantum numbers. Thus, we are looking for a decomposition of the form 

J'\^. y) = T.{y"Y'" f;l.,M^li'U^) , (le) 

The angular part of Z^" is given by the irreducible tensors of the four-dimensional rota- 
tional group Tj^-J,^^, which are orthogonal 

and satisfy the conditions: 

(n ) = T^^ in ) T^^ ) = 

Ti^LMyKt^^M = liCf\nyny^ (18) 

The measure duj in Eq. (|T7| ) relates to integration over the angles of unit vector n^, and 
Cf^ in Eq. ( [T^ ) are the Gegenbauer's (ultraspherical) polynomials. The polynomials 
Lnt{u) obey the condition: 

oo 

jdup,{u)L^,{u)K,,{u)=5^^,. 



The weight function pi{u) arising from the exponential term in ([l5| ) looks like 

ptiu) = M^e"", 

hence Lni{u) are the generalized Laguerre's polynomials. 
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The details of calculation of the currents J'^" 



bJln 



in Eq. (|T6D can be found in pa- 



per pi. As a result, the interaction term L2 takes the form 



^2 



^ bJln 

Cjg' 



G 



Jin 



A 



J, 



bJin 
fii...tMe 



[x) 



n\{e + n)\' 

V {xY 



(19) 
(20) 



" 'I 

'v 



2 



A 



{xY 



A2 



I 111. ..lie I 



IV (x) 

i 



A 



I}}. 



(21) 
(22) 



The doubled brackets in Eq. (^Tj) mean that the covariant derivatives commute inside these 
brackets. Form-factors Fne{s) are entire analytical functions in the complex s-plane, which 
is a manifestation of the gluon confinement. 

The classification of the currents will be complete if we will decompose J^a jf^ t^i "^ith 
J = V^A and £ > into a sum of orthogonal currents ^ail^^.fie with the different total 
angular momentum j = £ — l,i,£ + 1. Index a relates to the matrices = 7^ and 
F^ = 757a in Eq. (|2l|). This decomposition can be arranged by the following division 



bJin 
a,^\...fie 



i+1 

E-j-bJlnj 

j=l-l 



(23) 



where 



-j-bjinj 
-^a^l...^^l 



(£+1)2 ' atj.i...iii 



X q-bjln 
'JafJ.l'J p,pfj.2...fii 



e+1 



—V 7 



bJin 
bJin 



' a,fj.i...fii 



^5 7 



jbJin 



J 



bjen 

p, pp.2. ..fig 



3 = * 



(24) 



Symbol Vapi...pi ( ^41) denotes a cychc permutation of the indices (a/ii.../i£). Let sj be 
defined as 

Sp = ss = 0, Sv = sa = 1, 

then, using the orthogonality of the currents with different j, we can rewrite interaction 
term L2 as 



E E 

ajln j=\£—sj\ 



a J In I 



X] 
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where we have introduced the notation 



jy^"^ = J^J^I for J = S,P, i>0. 
fj.i...fj,i fj.i...fii ' ' — 

This form is equivalent to Eq. (|I0|), but now the interaction between quarks is expressed 
in terms of the nonlocal quark currents, that are elementary currents of the system in the 
sense of the classification over quantum numbers. 

For large Euclidean momentum the vertices V°"^^^ behave as l/(j9^)^+^/^. Therefore, 
only the "bubble" diagrams, shown in Fig. 1, are divergent. These divergencies can be 
removed by the counter-terms of the form —2X{x)TtVS. 

To avoid an unnecessary complication of notation, it is convenient to introduce con- 
densed index A/" enumerating the currents with all different combinations of the quantum 
numbers a, J, i, n and j. The renormalized vacuum amplitude Z takes the form 

Z = J c/dvac J DqDqexp I^J J d^xd'^yq{x)S~'^{x,y\B)q{y) 

+ E ^2^1 1 d'x [lAx) - TrVj^Sf"^ . (25) 

2.5 Bosonization 

By means of the standard bosonization procedure [|T], |^ applied to Eq. ( P^D the amplitude 
Z can be represented in terms of the composite meson fields HH^ 

Z = N Jl[D^MewS^^JJd^xdS<!>M{x)[{n-M^)6{x-y) 

- hlrK{x-y)]<l>^{y) + I,^,m}, (26) 



lint = J d'^xi J d'^X2hj^hj^,^j^{xi) [Tj^j^,{xi, X2) - Sj^j^>Ilj^{xi - X2)] ^M'{x2) 
If f ™ 

m=3 ^ •' •' k=l 

J rfdvacTr {Vv-i(a;i)5'(xi,X2 | B)...Vx,^{xm)S{x^, Xi \ B)} . 
Meson masses Mj^ are defined by the equations 

W(-M^) = 0, (27) 



■Ml...Mrr, 



1 + -2 



A / 
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where Ilj^{—M^) is the diagonal part of the two-point function r_VA^', which corresponds 
to the diagram shown in Fig. 2. a. The fields {Af = {a, J,i,n,j}) with j > satisfy 
the on-shell condition 

which excludes all extra degrees of freedom of the field, so that the numbers ^ and j 
can be treated as the 0(3) orbital momentum and total momentum, respectively P]. 
The total momentum j plays the role of an observable spin of the state with a given 
Ar={6,J,£,n,j}. 
The constants 



V = i/VnV(-M2,) (28) 

play the role of the effective coupling constants of the meson-quark interaction. 

The quark masses m/, the scale A (strength of the background field) and the quark- 
gluon coupling constant g are the free parameters of the effective meson theory (p^)-(p8|). 

In the one-loop approximation, the interactions between mesons with given quantum 
numbers AT = {b, J, i, n,j} are described by the quark loops like the diagram in Fig. 2.b. 
The structure of diagrams in Fig. 2 is the same as in the standard NJL-model, but in 
our case the quarks propagate in the vacuum gluon field, and the meson-quark vertices 
are nonlocal. Due to this nonlocality the quark loops are ultraviolet finite. The whole 
diagram is averaged by integration over the measure (icXvac- 

Figure 3 illustrates the central idea of the method of induced nonlocal currents, that 
has been realized in this section. An effective four-quark interaction is represented as an 
infinite series of interactions between the nonlocal quark currents characterized by the 
complete set of quantum numbers {b, J, i, n,j}. The form of the currents is induced by a 
particular form of gluon propagator. This new representation of the four-quark interaction 
generates an expansion of any amplitude into the series of the partial amplitudes with a 
particular value of the quantum numbers. Each partial amplitude is ultraviolet finite at 
any order of expansion over degrees of the coupling constant g. The composite meson 
fields in Eq. (^) are nothing else but the "elementary" collective excitations, that are 
classified according to the complete set of quantum numbers of the relativistic two-quark 
system. 

It should be stressed, that the model (p6D satisfies all demands of the general theory of 
nonlocal interactions of quantum fields [|rT[], which means that Eq. defines a nonlocal, 
relativistic, unitary and ultraviolet finite quark model of meson-meson interactions. 

Now we would like to test, how this formalism works in the meson phenomenology. 

3 Meson spectrum and weak decay constants 

Let us rewrite Eq. ( P?] ) in more detailed form 

A^ + G],,Ilbjinj{-Mlj,^^; mf, mj,- A) = 0. (29) 
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The function Hbjenj in Eq.(^) is given by the diagonal part in the momentum represen- 
tation of the tensor 



nJf:^,...(a:-y;m;,m^,;A) = 

J da.,,Tr [v^^^'::^^{x)S{x,y \ B)V^\'!:^^ (y) S {y , x \ B) 



(30) 



Relation (|29|) is the master equation for meson masses. The function IT can be calculated 
using the representations (^) and (|^) for the quark propagator and ( |?I| ) for the vertices. 
The only point, that requires a comment, is an averaging over the space directions of 
the vacuum field. Actually we have to average the tensors /^j^, ffj.ufap and so on. The 
generating formula looks as 



where (p, 6 are the spherical angles, J^i, is an anti-symmetric tensor, J^^, is a dual tensor 
and ± corresponds to the self- and anti-self-dual vacuum field. In particular this general 
representation gives: 



3.1 Light pseudoscalar and vector mesons 

First of all, let us fit the free parameters of the model, taking the masses of vr, p, K and 
K* mesons as the basic quantities. Below, we will sometimes use a symbol of a given 
meson instead of the corresponding set of quantum numbers (for example, vr instead of 



If Mjr, Mp, Mk and Mj^* are taken to be equal to the experimental values, then the masses 
niu and of the u and s quarks as the functions of A are defined by Eqs. (|3^) , (^) . Using 
m„(A) and ms(A) in (P^, we find the value of A, that provides a simultaneous description 





(31) 




2n^(-M2; m„, m„; A) = t{p{-Ml- m„, m„; A), 
2n/c(-M|.;ms,m„; A) = n/c*(-M|^,; m^, A), 
2n^(-M^; m„, m„; A) = n/^(-M|; m^, m„; A), 
/ = -9AVn.(-M2; m„, m,; A). 



(32) 
(33) 
(34) 
(35) 
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of the strange and nonstrange mesons. An optimal value of the coupling constant g is 
calculated by means of Eq. (pS]) . By this way we arrive at the values: 



198.28 MeV, = 412.96 MeV, A = 319.46 MeV, g = 9.96. 



(36) 



Solution (|36D is unique. 

It is well-known, that there should be a special reason, which provides a small pion 
mass and splits the masses of pseudoscalar and vector mesons. Breaking of chiral symme- 
try due to the four-quark interaction and two independent coupling constants for pseu- 
doscalar and vector mesons {gp ^ gv instead of our parameter g) play the role of such 
reason in the local NJL- model. As has already been pointed out, an interaction of quark 
spin with the vacuum field leads to the singular behavior of the quark propagator in the 
massless limit and generates a non-zero quark condensate, which indicates breaking of 
the chiral symmetry by the vacuum gluon field. Now let us illustrate, that in our case the 
same spin-field interaction is responsible for small pion mass and for the mass-splitting 
between P- and V-mesons. 

Polarization function Ilj (£ = 0, n = 0, J = P, \^) can be represented in the form 



nj(-M2;m/,m/,;A) 
A2 



47r2 



Tr^, / dti j dt 



+ 



where 



M^ Fy\h,t2,si,S2) 

A2 $4(ti,t2,Si,S2) 

2v{l-AvHit2)Fy\si,S2 

mtl,t2,Si,S2) 
^l(ti,t2,Sl,S2) 



dsi 

nifrnf. 



dsn 



l-si 

l + Si 



1-S2 



S2 



F^'\s,,S2) 



A2 



exp 



(l-S?)(l-si)$i(ti,t2,Sl,S2) 



\2fA2 



^{ti,t2, Si, S2) 



"^'2(^1,^2, Si, S2) 



(37) 



(38) 



$1 = 2viti + t2)[si^l + S2Q + SiS2[l + 4:vHM^i - 6 
$2 = 2f (ti + t2)(l + S1S2) + (1 + AvHit2)iSi + S2), 
FP = (1 + SiS2)[AiA2 + Av^ti - t2f^l^2SlS2], 



= 3 - siS2)AiA2 + 4t;^(ti - t2)"66siS2(l - 3siS2) 

Ai = [1 - AvHM^i - 6)]si + 2v{h + t2)6, 
A2=[l + 4:vHit2i^i - 6)]S2 + 2vih + ^2)6, 



1 



2 



{l + SlS2f 
2(l + SiS2) 



2 



1 _ 2 2 
1 - S1S2. 



(39) 
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Equations (p7D-(p^) shows, that singularity (1 — S] 



S2) ^, arising from the spin- 



field interaction (see the second line of Eq. (|D), is accumulated in the term of Eq. (^) 
proportional to the quark masses. Other terms are free from this singularity, although 
the spin-field interaction contribute to them either. This is due to the structure of the 
trace of the Dirac matrices. 

Let us compare a behavior of pion and p-meson polarization functions in the limit 

ruf = mfi = rriu <^ A. 



Using the singularity of the integrand in Eq. ( |37D at Si — > 1 and S2 — > 1, one can check, 
that the pion polarization function is singular in this limit and behaves as 



-M^;m„,m„; A) 



Tr, 



4f2A^ /-i/-! dtidt2 



exp 



r m2 1 



TT^ml Jo Jo $1(^1,^2,1,1) 

On the contrary, the p-meson polarization is regular at = and looks like 

np(-M^;m„,m„; A) = 
A2 



(40) 



An- 



Tr^^ dti dt2 dsi '^■52 jexp |^^^$(ti, ^2, Si, S2)| • (41) 



M2 )(ti, t2, si, S2) ^ 2y{l - Av%h)F^'>{s,, s^) 



A2 mti,t2,Si,S2) 

exp|^'^'(ii,i2,si,l) 



'^'2(^1,^2,51,52 
2v 1 



+ 



$2(^1,^2,^1,1), 



0(m„) 



and leads 



This difference appears owing to the form of factors -^2^'' and -^2^'' in Eq. ( 
to the inequality 

|n^(-M^;m„,m„;A)| > |np(-M2; m„, m„; A)|. 

This relation shows that the masses of pion and p-meson, satisfying Eq. (^), are strongly 
splitted and ^ when the quark mass goes to zero. A similar picture takes place 
for K and K* mesons, but since the strange quark mass is not so small, the effect is more 
smooth. 

Above consideration is illustrated in Table 0. The pion mass is much larger, and 
difference in the masses of pseudoscalar and vector mesons is smaller, if the spin-field 
interaction in the quark propagator is eliminated (compare the first and the last lines in 
the table). 

Thus, we can conclude, that the large splitting between the masses of P- and V- 
mesons is explained in our case by the spin-field interaction. This splitting is the reason, 
why Eqs. (^) have appropriate solution (^). 
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It should be noted, that the scalar polarization function Us differs from the pseu- 
doscalar lip only by the sign before mfmf in Eq. (0). Due to the above-mentioned 
singularity, the term proportional to mfirifi is leading, and is positive for a wide range 
of parameter values. As a result, Eq. (^) has no real solutions for the case of scalar 
mesons. It looks interesting, that the scalar {qq) bound states do not appear due to the 
same spin-field interaction, that diminishes the pion mass and provides nonzero quark 
condensate. 

Consideration of the SUpi'i) singlet and the eighth octet states shows an ideal mixing 
both for vector and pseudoscalar mesons. The masses of uo and 0, calculated with the 
parameter values (^), are 

M^ = Mp = 770 MeV, = 1034 MeV, 

which is in a good agreement with the experimental values. An ideal mixing of the 
pseudoscalar states is not the case, that can provide an appropriate description of rj and 
rj' mesons. It is well known, that the problem of f] and 1]' masses can be solved by taking 
into account another Euclidean gluon configuration - the instanton vacuum field [0 . The 
instantons can be incorporated into our formalism without any principal problems, and 
we hope to realize this idea in forthcoming publications. 

Now let us consider the weak decays of vr and K mesons. In the lowest approximation, 
an amplitude of the decay P W is given by the diagram in Fig 2.c. The weak decay 
constant fp is defined by the standard formulas 

Ap^Uk, k') = t^}ChpFie)^p{k)kJ{k'){l - ^,h^u{k + k'), 

fp = hpF{-M^), (42) 

where the meson-quark coupling constant hp is calculated via Eq. (|28D, and /C is the KKM 
matrix element corresponding to a given meson. For arbitrary pseudoscalar meson the 
diagram in Fig. 2.c gives the following expression for fp 

r _ 1 f f dtdsids2{l + S1S2) /I - sA 3^ /I - S2\ 

^4^ ^"7770 [2vt{l + S1S2) + si + 82]^ U + sJ \l + S2 



si + 2vt[l - + sj)] S2 + 2vt[l - ^2(1 + si 
\ -ji 



I -si ■ - I- si 



^ ^ S1S2 + 2vt{siil + S2il) /43X 

2vt{l + S1S2) + Si + S2 

The singularity of the integrand of Eq. ( |43D at si — > 1 and S2 ^ 1 appears from above- 
mentioned spin-field interaction in the quark propagator and plays the main role in regu- 
lating the value of fp for the light mesons. Calculation of pion and kaon decay constants 
by means of Eq. (^) with the values of parameters (|36D gives 

= 126 MeV, fK = 145 MeV. 
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Note, that the couphng constants h.,^ and hx depend on the meson mass, quark masses 
and parameter A (see Eq. (|28| ) and Table |). 

One could get a definite impression, that simultaneous description of the masses of vr, 
K, p, K*, uj and mesons, and quite accurate values of f.^ and fx are obtained mostly 
due to the breakdown of chiral symmetry by the spin- field interaction (see also Eq. (0)). 

In order to clarify a status of this impression, one needs to investigate the chiral limit 
of the model. Whether or not, and in what form the Goldberger-Treiman and Oakes- 
Reiner-Treiman relations are fulfilled in this limit? Naively, it seems that the chiral limit 
corresponds to the case ^ A. However, the situation is much more complicated. 

Just for illustration of this statement, consider the pion mass in the limit ^ A. 
Integrals in Eq. (1^) can be evaluated, and we get the following asymptotic form of 
equation (E^) for pion mass 



16A6 



9 



exp 



Ml 
8vA^ 



Ml 



n 2 



exp 



0. 



(44) 



For any values of g there is a real positive or negative solution to Eq. (^^. The pion 
mass is equal to zero, if the values of m„, g and A satisfy the relation 



m 

X2 



2g_ 

IStt 



For a fixed value of g and m^/A 0, the solution to Eq. (^if ) is purely imaginary 
and behaves as 

In ln(A/m^j) 



A 3 V nT'ti 



y^ln(A/m„) 



which has nothing physically reasonable interpretation, but just indicates, that the limit 
g = const, and rriu/A — is ill-defined. 

/^From our point of view, a correct transition to the chiral limit has to be based on a 
simultaneous changing of rriu, g and A as functions of an actual physical parameter like 
the temperature or particle density. The quark masses, coupling constant and the vacuum 
field strength as the functions of this parameter have to be extracted from consideration 
of QCD dynamics at nonzero temperature and density. Unfortunately, this is a quite 
complicate problem, and we leave it for further investigations. 

The main result of this subsection is very simple: we have demonstrated by the explicit 
model calculations, that the spin-field interaction contained in the quark propagator (in 
presence of the homogeneous (anti-)self-dual vacuum gluon field) can be responsible for 
observable masses of the light pseudoscalar and vector mesons (by exception of r], rj'), and 
for the values of the weak decay constants of pion and kaon. All numerical results are 
given in Tables |l| and 0. 
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3.2 Regge trajectories 

It has been shown in our previous paper P, that the spectrum of radial and orbital 
excitations of the light mesons is asymptotically equidistant: 

Mlj.n = ^ In (^) • ■ n + O(lnn) , for n » £, (45) 

MIj^^ = ^ In 5 ■ A^ ■ ^ + 0(ln£), for £ > n. (46) 

Technically this result is based on the exponential behavior of the quark propagator (^) 
and vertex function F^i p^) in the Minkowski region (p^/A^ — > — oo) like 

HfiP I 5) - O (-P {J^}) , (f) - O (exp ||^}) , (47) 

and on the specific dependence of the coupling constant G jgn (see Eq. (|l^)) on the orbital 
and radial quantum numbers £, n, arising from the decomposition of the bilocal quark 
currents over the generalized Lagguerre polynomials, which is determined in its order by 
the form of gluon propagator (|^) (for details see sect. ^.41 ). In general words, the Regge 
character of the spectrum is determined in our model by the confining properties of the 
vacuum field. 

Numerical calculation of masses of the first orbital excitations of vr, K, p and K* 
mesons by means of Eq. (|29|) with the parameters (36) gives the masses shown in Table |[ 



The super-fine structure of the excited states of p and K* mesons coming from classi- 
fication of currents over total momentum (Eq. (|23|) ) is qualitatively correct. Super- fine 
splitting of the levels with £ = 1 is not very large. 

3.3 Heavy quarkonia 

Exponential behavior of the quark propagator and vertices (0) is responsible for the 
following relation between the masses of heavy quarkonia Mqq and heavy quark mq in 
the leading approximation 0: 

Mqq = 2mQ, for uiq > A. 

Now let us calculate the next-to-leading term in the mass formula. In other words, we 
have to solve Eq. ( |29D with the polarization function IIj defined by Eq. (^) with 

= = "^Q > A, Mqq = IVflQ - Aqq (48) 

in the next-to-leading approximation over I/itiq. Since the masses of quarks are equal to 



each other, we have (see Eq. (|T3D) 

^1 = 6 = 1/2, 
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which means that the composite quarkonium field ^qq{x) is locahzed at the center of 
masses of two heavy quarks (in the Euchdean four-dimensional space). It is convenient 
to transform variables Si and S2 in Eq. (P7|): 



ri = (si + S2)/V2, r2 = (si - S2)/V2. (49) 

The term with ^3'^^ in Eq. ( |37D does not contribute to the leading and next-to-leading 
behavior of the integral and can be omitted. After the transformation we arrive at the 
expression: 

Ilj{-M^;mQ,mQ;A) = 

rriQ r ( ri/v^ f^i rV^-ri \ 

-—^Tr^ // dtidt2 / dri \ dr2+ / dri / dr2\ ■ 

47r^ JJO Wo J-ri J1/V2 J~V2+ri I 



Ri'\n,r2) 



+ 



R[''\ti,t2,ri,r2) 

V2{ti,t2,ri,r2) ' (2 - (ri - r2)2)(l - (n + r2)2)v9|(ti,t2,ri,r2) 

O 



ew{^vih,t2,ri,r2) 



where 



^ y^i(ti,t2,ri,r2) _ ^ (V2 + ri)^ -r| 

= V2v{ti + t2)ri + ri - rj 
^2 = v{h + t2){2 + r\ - r2^) + 72(1 + ^vH^t2)r\, 



(50) 



(51) 



-(2 + rl - tI)[A,A2 + - t2f{Tl - vDl 



R'-P = Y^[(6 - + rl)A,A2 + i;^(ti - t2f{rl - rl){l - 3r? + Sr^)], 

^1 = ri - r2 + V2t;(ti + t2), A2 = n + r2 + V2t;(ti + ^2), 

l + rl) 



J^2 



{2 + rl-rl)\ R^r=A-{r',+rl^' 



An asymptotic value of the integral over r2 in the limit M ^ A can be evaluated by the 
Laplace method. One can check, that the function has a maximum at the point r2 = 
for any values of ri, t\, t2: 



_d_ 

dr2 

dr2 



V9(tl,t2,ri,r2)|r2=0 = 0, 



8v^ri Av{h + t2) + ^ri{l-v^{ti-t2f) 



M2 (2 - ) 



2\2 



-^^{ti,t2,ri,r2)\r^=o - 

ri[ri + 2^v{ti+t2)] 2m^^^V2 + ri 



V2iti,t2,ri,Qi) 



V?(ti,t2,n,0) 



V52(ii,i2,n,o) 



M2 72- 



ri 



<0, 
(52) 
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which means that the leading terms can be obtained by evaluating the Gaussian integral 
over r2. Furthermore, one can see, that the largest value of the function ip{ti,t2,ri,0) in 
the interval ri G [0, \/2] corresponds to ri = for any ti, t2, moreover: 

— ^(ti,t2,ri,0)|.,=o = -^(^^-lj <0, (53) 

and, in the leading approximation, the integrand is reduced to the exponential function 
in ri. Using Eqs. (^) and (p^D, and taking into account Eq. (^), one can integrate over 
r2 and ri with the result: 

2 3A^ ff^ dtidt2 



UA-M'; ,n„ ,n,; A)| = - g + 0{^]. (54) 




Integrating over ti and t2 in Eq. ( p^ and substituting the result to Eq. (pQ]), one can find 

A TTy/n \T^Q/ 

where Cp = 1/9, Cy = 1/18 (see Eq. (|I9|)). It should be stressed, that the difference in 
the constants 

originates from the Fierz transformation of the Dirac matrices in the interaction term L2 
in representation (0). Relation (p6D means that the vector quarkonium state is always 
heavier than the pseudoscalar one. 

The results of numerical calculation of the masses of different heavy quarkonia states 
are summarized in Tables ^ and |^. The parameters A and g are equal to the values (^) 
fitting the light meson masses, and = 1650 MeV, rrih = 4840 MeV. An agreement with 
the experimental values is rather satisfactory. The super-fine splitting {xco, Xci, Xc2 and so 
on) is very small, since it is regulated by the terms ©(l/mg) in Eq. (pOj) . Its description 
is qualitatively correct. The splitting is generated in our model by dividing the quark 
currents with 7^ and i > into the antisymmetric, symmetric traceless and diagonal 
parts (see Eq. (|23|) ), which extracts the states with different total angular momentum, 
mixed in the currents gMjaT^^^ ^^Fniq. 

We conclude that the correct description of the heavy quarkonia in our model is 
provided by the specific form of nonlocality of the quark and gluon propagators induced 
by the vacuum field, localization of meson field at the center of masses of constituent 
quarks and by a separation of the nonlocal currents with different total momentum. In 
general, the spectrum is driven by the rigid asymptotic formulas (ESI) and 
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3.4 Heavy-light mesons 

Another interesting sector of meson spectrum is the heavy-hght mesons, characterized by 
a rich physics [|l^, 0|. In this subsection we will consider the masses and weak decay 
constants of heavy-light mesons. First of all, let us get the asymptotic formulas in the limit 
of infinitely heavy quark. Namely, we have to investigate the behavior of the polarization 
function nj(— M; mg, m^; A) (Eq. and the weak decay constant fp (Eq. in the 
case 

TTLf = TUq ^ a, = ~ C'(A), 

if = "1^ = 1 + OimJrriQ), = = OimJniQ). (57) 

' niQ + ITLq ^ ^ ' niQ + mq 

Equations (|57D indicate, that in the heavy quark limit the composite meson field ^Qq{x) 
is localized at the point, in which the heavy quark Q is situated. 

Let us show, that in the limit (^7|) the leading and next-to-leading terms of the solution 
to Eq. (|2|) read 

Mq^ = mQ + aJ/) + 0(l/mQ), (58) 

where the next-to-leading term Aq^ does not depend on the heavy quark mass mq. This 
term is a function of a light quark mass niq and coupling constant Gjoo (see Eq. (p!9|)). 

Omitting the term with which does not contribute to the leading and next- 

to-leading behavior of the integral and taking into account conditions (^), one can 
rewrite Eq.(p7D in the form 

nj(-M^;mQ,mg; A) = ~J:^^^^JJJJq dtidt2dsids2 ^ _^ ^ 
;i - Av%t2)Y{tu S2)T['\si, S2)s,M^ , 



[s,X{ti,t2,S2)+Y{ti,t2,S2)]^ 

Ty\si,S2)mQmq 

:i - Sl){l - sl)[siX{tu t2, S2) + Y(ti, t2, S2)Y 



exp|^^0(ti,t2,Si,S2) !> + ( ) , (59) 




where 



siY{ti,t2,S2) ruQ X + si 



SiX{h,t2,S2)+Y{h,t2,S2) 2M2 1-Si 

X = l + 4v%t2 + 2v{ti+t2)s2 

Y = 2v{ti + t2) + (1 + 4v\t2)s2 
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n 



1 + S1S2, 

(l+SlS2)^ T 



(V) 



-(3 - S1S2) 

1 22 
1 - 



One can check, that for any ti, t2 and S2 the function 0(^1,^2,-31,32) has a maximum at 



rrir 



A' 



0(^l,^2,s^^s2) = 

^0(tl,t2,Sl,S2)U=s- 

Therefore, we can write 

nj(-M2;mQ,mg;A) = 
1 



M2 

^^(^1,^2, S2 



2 \ 2 



2X(ti,t2,S2) 



M2 



o 



2X(ti,t2,g2) 
>'(tl,t2,S2) 



2 \ 2 



(60) 



47r2 



Tr, 



dtidt2ds2 



I-S2 



4iiA2 



M2 



\ ^q\' y(M2,g2) 

4t;A2 \ My X(ti,t2,S2) 



1 + S2 

1 - 4:vHit2)s'^^''M^ 

Y%h,t2,S2) 

dsi exp 



exp 



+ 



rnqmy 



sl)Y^{h,t2,S2) 



2vA^Y{h,t2,S2) \ 



O 



■ 



(61) 



Integrating out the variable Si in Eq. (0), substituting M = ttlq + A^^^ to the resulting 
expression and using Eq. (^9]), we arrive at the equation 



1 = g'Cj 



^^qI X{h,t2,S2) m. 



1 



Tr, 



(27r)3/2 



A 



^dtidt2dS2 [il-S2)/{l + S2)]~- 

[X{t,,t2,S2)Y{t,,t2,S2)f' 
2 

^-(^1^2,52) 



exp 



a(J) 



VA^ X{t,,t2,S2) 



+ 



-) 



(62) 



Equation (^) describes dependence of Aq^ in mass formula (|58|) on the coupling constant 
g, the light quark mass and vacuum field strength B (A). There is a single real solution 
to Eq. (|6^) for any positive g, rriq and A. In particular, for the values (^) we get 



Ag;) = 20 MeV, 



A 



iP) 



63 MeV, 



a(v) . 



155 MeV, 
: 191 MeV. 
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As is seen from Eq. (|62D, the difference between pseudoscalar Ag^^ and vector Aq^"* is due 
to the constant Cj, that appears from the Fierz transformation of the Dirac matrices. 
This is the same situation as in the case of heavy quarkonia (see Eq. (p5|)). 

Table ^ demonstrates reasonably good agreement between the experimental data and 
the masses of the heavy-light mesons calculated by means of Eq. (|29| ) with the parame- 
ters (|36|) . The masses of heavy quarks are the same as in the description of the heavy 
quarkonia (see Table |l|). 

Now let us turn to the calculation of the weak decay constant for the pseudoscalar 
heavy-light mesons. Under conditions (^), the integral over si in Eq. (|43|) can be evalu- 
ated by the Laplace method. The result is 



ff 




(63) 



where 



Af = Tr,, 



A 



(P) 

Qq 



(27r)3/2 



f f^dtdS2[{l-S2)/{l + S2)]^ 



[S2 + 2vtf^ [1 + 2vtS2] 



3/2 



A 



+ 



TUq 1 



2vtS2 



A 



s2 
"^2 



exp 



A 



Qq 



S2 + 2vt 



VA^ 1 + 2vtS2 




(64) 



(p) 

and the difference Aq^ between the masses of heavy quark and heavy-light meson is 
given by Eq. (0). The procedure for obtaining Eqs. (^) and (p^ ) is very similar to the 
calculations providing Eq. (|62|). To get the final formula for fp, the asymptotic form 
of meson-quark coupling constant hp has to be defined in the case (0). Performing 
calculations analogous to that ones, which lead to Eq. (|6^), we arrive at 



A: 



a{P) 
^Qq 




^ 2(27r)3/2A 

A^ ^ X{ti,t2,S2) ruq 



^dtidt2dS2 [(l-g2)/(l + S2)]^- 
[X(ti,t2,S2)F(tl,t2,S2^l'/' 



A 



A 



exp 



a(J) 
^Qq 



Y{ti,t2, S2 



VA^ X{ti,t2,S2 




— ■ (65) 



One can see, that Eqs. (|63|)- (|65|) give the following asymptotic relation in the heavy quark 
limit (13) 

fp = — (66) 



23 



where Aj and Ah do not depend on the heavy quark mass in the leading approximation 
over A/rriQ, as is indicated in (|64D , (|65D . Relation (|66D agrees with the accepted notion 
about behavior of the weak decay constants of the heavy-light mesons |T^. Results of 
numerical calculation of the weak decay constants for different pseudoscalar mesons are 
given in Table ^ 



4 Discussion 

For conclusion, we would like to point out several problems, that require more profound 
studying. 

We have assumed from the very beginning, that the nonperturbative QCD vacuum 
is characterized by a nonzero background (anti-) self-dual homogeneous field. In other 
words, a minimum of the QCD effective potential (the free energy density) for this gluon 
configuration is assumed to be at nonzero field strength. Different estimations of the 
effective potential indicates, that this situation can be realized Although these 



estimations cannot be used as a basis for more or less rigorous proof, they underline the 
key role of the asymptotic freedom in forming the effective potential for an homogeneous 
gluon field. Just the asymptotic freedom is the distinctive feature of nonabelian gauge 
theories like QCD. Gluon self-interaction is manifested also in the nontrivial form of gluon 
propagator (^). Although the background field is quasi-abelian, the nonabelian nature of 
the gluon field plays the crucial role for above considered model. 

In order to clarify the basic assumption of this paper, one needs to get a reliable 
nonperturbative estimation of the free energy density or effective potential of QCD for the 
background field under consideration. Lattice calculations seem to be the most promising 
approach to this problem. 

In this paper, we have demonstrated how the singularity l/mf of the quark propagator 
affects the masses and weak decay constants of light mesons. However, more detailed 
consideration of the chiral symmetry breaking by the background field is needed. This 
can be achieved by investigating the Dirac equation in presence of the homogeneous (anti- 
)self-dual field. Another source of violation of the chiral symmetry is the effective four- 
quark interaction. The divergent diagram in Fig. 1 should play the key role in studying 
of this mechanism of the symmetry breaking. One can expect that additional breakdown 
of chiral symmetry by the four-quark interaction could diminish the quark masses. They 
come into our formalism as the current masses, but their fitted values (see Table |1|) are 
close to the constituent quark masses rather than to the current ones. 

One could see, that the coupling constant g in Table |1| is rather large. A possible 
origin of this unpleasant feature could be covered in elimination of some terms of the 
gluon propagator (|^) (for more details see In other words, some truncations in the 
gluon propagator was compensated by the rising of the coupling constant. This point also 
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has to be investigated carefully. 
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Table 1: Parameters of the model. 



niu (MeV) 


rud (MeV) 


m, (MeV) 


rric (MeV) 


nib (MeV) 


A (MeV) g 


198.3 


198.3 


413 


1650 


4840 


319.5 9.96 



Table 2: The masses(MeV), weak decay constants (MeV) and meson-quark coupling 
constants h of the light mesons. M* - calculation without taking into account the spin- 
field interaction. 



Alt^sou TT A' A* L^' o 

M 140 770 496 890 770 TosI 

M'^^P 140 770 496 890 786 1020 

fp 126 - 145 - 

/^^P 132 - 157 - 

h 6.51 4.16 7.25 4.48 4.16 4.94 

M* 630 864 743 970 864 1087 
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Table 3: The masses (MeV) of orbital excitations of tt, p and K* mesons. Super-fine 
structure of the I — 1 excitation of p and K* is shown (£ is the orbital momentum and j 
is the total momentum (an observable spin) of a state). 



Meson 


e 


J 


M 


^exp 


TT 








140 


140 


h 


1 


1 


1252 


1235 


K 








496 


496 


/s:i(1270) 


1 


1 


1263 


1270 


p 





1 


770 


770 




1 





1238 




ai 


1 


1 


1311 


1260 




1 


2 


1364 


1320 


K* 





1 


890 


890 




1 





1274 




ii'i(1400) 


1 


1 


1342 


1400 


K* 


1 


2 


1388 


1430 



Table 4: The spectrum of charmonium. 



Meson ?7c J/tp Xcp Xc2 V^' j'" 

n 1 2 

£ 1110 

j 0101211 
M (MeV) 3000 3161 3452 3529 3531 3817 4120 
M^^P (MeV) 2980 3096 3415 3510 3556 3770 4040 



Table 5: The spectrum of bottomonium. 



Meson T Xb„ Xh Xh-, ^' Xb,, Xb, ^" 

n 1 1 1 1 2 

e 011101110 

j 101210121 

M (MeV) 9490 9767 9780 9780 10052 10212 10215 10215 10292 
M^^P (MeV) 9460 9860 9892 9913 10230 10235 10255 10269 10355 
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Table 6: The masses and weak decay constants (MeV) of heavy-light mesons. 



Meson 


D 


D* 


Ds 


d: 


B 


B* 


Bs 


b: 


M 


1766 


1991 


1910 


2142 


4965 


5143 


5092 


5292 




1869 


2010 


1969 


2110 


5278 


5324 


5375 


5422 


fp 


149 




177 




123 




150 





A 

FIG.l. Divergent bubble diagram. 
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a 



FIG. 2. Diagrams describing different processes in effective nonlocal meson theory in 
the lowest (one-loop) order. 
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aJini 



■laJinj 



aJinj 



Gluon propagator in the vacuum field 




Local quark current 



> 



Nonlocal quark currents 



FIG. 3. The decomposition of the four-quark interaction. 
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